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ABSTRACT 

-  V.  / 

Wtrconsider  the  estimation  of  a  collection  of  location 
parameters  when  it  is  believed,  a  priori ,  that  their  ordering 
is  known.  The  least  squares  and  least  absolute  deviations 
estimates  subject  to  this  ordering  restriction  have  been 
studied  in  the  literature.  We  seek  robust  estimators  which 
perform  well  for  a  broad  range  of  distributions.  The  results 
of  a  Monte  Carlo  study  and  a  study  of  computation  algorithms 
are  discussed. 


(^This  research  was  sponsored  by  the  Office  of  Naval  Research 
under  ONR  contract  N00014 - 80 - C- 032 2  . 
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Parts  of  this  work  are  taken  from  this  author's  doctoral 


dissertation  written  at  the  University  of  Missouri -'Ro  1  la. 
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1.  INTRODUCTION  AND  SUMMARY.  We  consider  the  estimation  of 
k  location  parameters  when  it  is  believed  that  they  arc 
nondecreasing.  Let  0^  <_  62  £  .  ..  £  6^  denote  the  parameters 
and  suppose  that  independent  random  samples,  X^  ,  j  =  1  ,  2  ,  .  .  .  , 
n^,  i=l,2,...,k,  are  available.  Brunk  (1955)  obtained  the 
maximum  likelihood  estimates  of  nondecreasing  normal  means, 
which,  of  course,  minimise 

^i=l^j=l(Xij'6i)2  subject  to  6j  <  62  <  ...  <  6k> 

There  are  several  algorithms  for  computing  these  estimates, 
but  we  emphasize  the  pool  adjacent  violators  algorithm  ( PAVA ) . 
(For  a  detailed  discussion  of  such  algorithms,  see  Section  2.3 
of  Barlow  et  al.  (1972).)  If  the  sample  means,  X-,  are 
nondecreasing,  then  they  are  the  restricted  least  squares 
estimates.  If  not  there  is  a  violation,  that  is  >  Xj+J 
for  some  i,  then  the  corresponding  samples  are  pooled  and  the 
two  sample  means  are  replaced  by  the  mean  of  the  pooled 
sample.  Next  the  resulting  k-1  sample  means  are  considered 
with  the  understanding  that  once  two  samples  are  pooled  they 
must  remain  together.  This  process  is  continued  until  a 
nondecreasing  set  of  means  is  obtained.  If  the  ith  sample 
has  not  been  pooled  with  another,  then  the  estimates  of  0. 
is  X^,  but  if  the  ith  sample  has  been  pooled,  then  the  estimate 
is  the  mean  of  the  final  pooled  sample  containing  the  ith 
sample . 

As  would  be  expected,  these  restricted  least  squares 
estimates  are  unduly  affected  by  extreme  observations  and  so 
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robust  estimators  would  be  desirable  in  many  situations. 

Robertson  and  Waltman  (1968)  derived  the  least  absolute 

deviations  estimates,  in  particular  they  obtained  the 

values  which  minimize 
.  n . 

z-i  =  l£j  =  ll  ij  l1  J  1  -  2  -  -  k 

They  also  showed  that  if  one  uses  the  sample  medians  as  the 
initial  estimates  (in  place  of  X • )  and  computes  the  median 
of  the  pooled  samples  rather  than  the  mean,  then  the  PAVA 
can  be  used  to  compute  these  estimates  also.  It  is  well- 
known  that  the  sample  median  is  not  very  efficient  for  a 
normal  model,  but  this  does  suggest  an  approach  for  obtaining 
robust  estimates  of  the  0^.  Using  the  notation  from  Huber 
(1981),  let  T  be  a  location  statistic,  such  as  the  mean, 
median  or  an  M-estimate,  which  may  be  thought  of  as  a 
functional  defined  for  all  empirical  distribution  functions 
(EDFs).  One  could  obtain  nondecreasing  estimates  of  the  0. 
by  employing  the  PAVA  and  using  T  to  obtain  initial  estimates 
as  well  as  the  estimates  from  pooled  samples.  We  denote 
such  estimates  by  _<  0^  £  .  .  .  £  0^.  (While  these  estimates 
depend  on  T,  we  do  not  indicate  this  in  the  notation.) 

Robertson  and  Wright  (1980)  considered  computation 
algorithms  for  order  restricted  estimates  and  found  that  it 
is  desirable  for  the  estimator,  T,  to  have  the  Cauchy  mean 
value  property  (CMVP) .  T  has  the  CMVP  provided  its  value 
for  a  pooled  sample  is  between  the  T  values  for  the  two  samples. 
That  is,  if  P  and  F  are  EDFs  based  on  samples  of  size  n  and 
m  respectively,  then 
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TfFn)AT(Fm)  -  T((nFn+mFm)/(n+m))  <  UF^VTO-j  . 

An  estimator  with  this  property  will  be  called  a  Cauchy 
mean  (CM) .  If  T  is  a  CM,  then  the  PAVA  yields  the  same 
estimate  independent  of  the  order  in  which  violators  are 
pooled  and  these  estimates  agree  with  those  obtained  from 
the  inax-min  formulas  which  are  discussed  there.  Robertson 
and  Wright  (1974)  have  shown  that  if  T  is  a  consistent  CM, 
then  the  estimates  obtained  from  the  PAVA  also  are  provided 

n.  00  for  each  i. 

1 

In  Section  2,  several  classes  of  robust  estimators  are 
studied  to  determine  if  they  are  CMs.  Since  so  few  of  those 
considered  in  the  literature  are,  the  use  of  non -CMs  is 
discussed  further.  Section  3  describes  the  results  of  a 
Monte  Carlo  study  of  robust  estimators  for  ordered  location 
parameters.  If  one  wants  estimators  that  perform  well  over 
the  range  of  normal  to  Cauchy  errors,  then  Gastwirth's 
estimator  or  the  trimmed  mean  which  trims  25%  on  each  side 
are  recommended.  The  trimmed  mean  is  a  little  more  efficient 
for  the  normal  model  and  Gastwirth’s  estimator  is  more 
efficient  for  the  Cauchy.  The  latter  is  also  easier  to  compute 
If  very  heavy  tailed  distributions  are  not  a  concern,  we 
recommend  the  Huber  with  c  =  1.5  (c  is  defined  in  Section  2) 
which  has  the  advantages  of  the  CM$ . 

2^  ROBUST  CMs  AN 0  COM PUT ATJ ON  ALGOR l TUNIS .  In  this  section, 

we  consider  the  types  of  location  estimators  discussed  in 
Andrews  ct  al.  (1972)  to  determine  which  of  those  would  he 
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appropriate  for  computing  order  restricted  estimates. 

An  M-estimate  is  a  solution,  =  T(F  J,  to  an  equation 
of  the  form 

Ij  =  l^( (xj  ~Tn)/s)  =  °» 

where  <p  is  an  odd  function  and  s  is  estimated  independently. 

(One  can  also  estimate  s  simultaneously,  but,  as  we  shall 
see,  this  may  yield  an  estimator  which  does  not  have  the 
CMVP.)  Hampel  and  Andrews  have  proposed  some  redescending 
ip  functions  (cf.  Andrews  et  al.  (1  972)),  but  Hogg  (1  979) 
points  out  that  the  M-estimates  corresponding  to  these 
tp  functions,  as  well  as  Tukey's  biweight,  may  possess 
convergence  problems  when  solving  iteratively.  Since  such 
an  iterative  procedure  must  be  implemented  several  times  when 
computing  order  restricted  estimates,  these  <1  functions  were 
not  considered  further.  Huber  (1964)  proposed 

^(x)  =  x  for  |  x  J  _<  c  and  iJj(x)  =  c  sgn(x)  for  |x|  ■  c,  (2) 

with  c  a  fixed  positive  constant.  A  common  choice  for  s  is 
median(|x^-median(x^) | ) / . 67  54  , 

however,  it  is  not  difficult  to  construct  examples  to  show 
that  if  s  is  computed  this  way,  varying  with  the  sample, 
then  the  resulting  estimator  need  not  be  a  CM.  (An  example 
is  given  in  Magel  (1982).)  In  many  situations,  the  k 
populations  are  assumed  to  differ  only  in  location  and  so 
one  could  use  a  fixed  value  of  s,  namely, 

s  -  med ian ( | x .  . -medi an (x . . :  j  =  1 ,  2  , . . .  , n  . ) | :  i  =  1  ,  2  , . . .  ,  k ) / . 67 
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It  is  easy  to  show  that  the  solution  set  for  this  ip  must  be 
a  nonempty,  closed  interval  (possibly  one  point)  and  so  we 
define  the  estimator  to  be  the  midpoint. 

Remark.  Huber '  s  M-est  imator  (ip  given  by  (2)J  with  a  fixed  s 
is  a  CM . 

Proof.  Without  loss  of  generality  we  assume  s  =  1.  Let 

T  =  T ( F  )  ,  T  =  T (F  )  and  T  T ( F  .  ,  where  F  and 

n  n  m  m  n+m  1  n+m) ’  n 

Fm  are  the  EDFs  corresponding  to  two  samples  and  F  m  = 

(nF^  +  mFm)/(n+m).  Let  the  solution  set  for  be  t a  j  » b  j  1  » 

for  Tm  be  [a7,b2]  and  for  Tn  +  m  be  [a^.b^].  We  also  assume 

that  the  samples  have  been  labeled  so  that  a^  £  a^-  Now 

the  left  hand  side  (lhs)  of  (1)  can  be  written  as 
’ 

n  (x-Tn)  dFn  (x)  and  we  see  that 
'l'(x-t)  dFn  (x)  a  jij>  (x  - 1 )  dFm(x)  1  |iKx-t)dFn+m(x) 

(3) 

<  tp(x-t)dF  (x)v  <Hx-t)dF  (x)  . 

—  j  n  m 

Setting  t  =  a^-£,  e  >  0,  the  middle  term  of  (3)  is  positive 
and  so  one  of  the  expressions  in  the  rhs  of  (3)  must  be 
positive.  Because  a^  £  a2  this  implies  that  a,  -  c  <  a9. 
Letting  e  -*  0,  yields  a^  £  a2>  Setting  t  =  a^  in  tJic  middle 
term  implies  that  one  of  the  terms  in  the  lhs  is  nonpositive 
or  a^  ^  a^.  So  if  T  were  chosen  to  be  the  lh  endpoint  it 
would  have  the  CMVP.  By  symmetry  considerations,  the  same 
can  be  shown  for  the  rh  endpoints.  If  b^  £  b^  ^  b2>  then 
(a^  +  b^)/2  £  (a^  +  b^)/2  _<  (a2  +  b2)/2.  If  b-,  <  b ,  «-  bj,  then 
[  a  2  ,  b  2  |  C  l  a  A  ,  b  1  ]  ,  in  which  case  |  ,  b  ^  ]  =  |a,,b.,J,  because 
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"«  t  "a  , 

L. 

’ 

(n  +  m)  ^(x-t)dF  +m(x)=n  (x  - 1 )  dF^  (x)  +m  ip  (x- t j  dFffl  ( x)  =  0  a?^t^b? 

-'()  h?-t. 


This  remark  applies  when  s  is  fixed,  but  in  the  above 

definition  of  Huber's  estimate  s  may  vary  with  the  k  samples. 

However,  for  a  fixed  set  of  variables  X..,  s  has  a  fixed 

ij 

value  and  if  that  value  is  used  to  compute  the  6.,  then  the 
order  in  which  violators  are  pooled  will  not  affect  0, 
the  max-min  formulaes  will  give  the  same  0.,  and  the  norm 
reducing  property  of  Robertson  and  Wright  (1974)  will  imply 
that 

max. „ .  .  I  6  .  -  6  .  I  <  max .  . .  ,  I T  -0.1. 
l^i <k 1  i  l1  —  l£i_<k  i 

So  if  the  error  distribution,  ie.  the  distribution  of 

X..-0.,  is  symmetric  about  zero,  then  0\  -+•  0 .  for  each  i 
ij  i  ’  ’ll 

provided  n^  •*  00  for  each  i. 

Another  large  class  of  location  estimators  arc  the  linear 
combinations  of  order  statistics,  or  L-estimators .  This 
class  includes  trimmed  means,  adaptive  trimmed  means, 
Gastwirth's  estimator  and  Tukey's  Trimean,  as  well  as  the 
mean  and  median.  Leurgans  (1981)  has  shown  that  there  are 
only  three  basic  types  of  L-estimators  which  possess  the  CMVP: 
the  mean;  weighted  midranges,  that  is  a  weighted  average  of 


the  smallest  and  largest  order  statistics,  wx ( ] j  +  (l-w)x(nj 

where  x , . .  <  x,_.<  ...  <  x,  .  are  the  order  statistics  from 
(1)  -  (2)  -  -  (n) 

a  sample  of  size  n  and  0  £  w  _<  1 ;  and  weighted  percentiles, 


x ( [ np ] + 1 )  if  np  is  not  an  inte«er’  and  wx(np)  + 


( 1  -  w )  x 


(np+  1 ) 


8 


if  np  is  an  integer.  Hence,  the  commonly  used  robust  I.  estimators 
are  not  CMs .  Hogg  (1967)  proposed  an  adaptive  estimator  which 
uses  various  L-estimators  depending  on  the  value  of  the 
sample  kurtosis.  Two  of  these,  the  trimmed  mean  and  the  outei 
mean  are  not  CMs. 

estimators  derived  from  rank  tests  for  a  shift  are  called 
R- est imators .  The  Hodges  -  Lehmann  estimator,  a  popular  member 
of  this  class,  is  defined  to  be  the  median  of  pairwise  averages, 
med ( (x ^+xj ) / 2 ) .  Magel  (1982)  gives  an  example  which  shows 

that  this  estimator  does  not  possess  the  CMVP  and  the  same 

■> 

example  shows  that  this  is  true  if  one  considers  all  n“  pairs, 
just  those  with  i  <.  j  >  or  just  those  with  i  <  i.  The  folded 
medians  comprise  a  closely  related  class  of  estimators  which 
depend  on  the  averages  of  symmetrically  placed  order  statistics, 
ie.  (xtl)+X(n))/2,  (x(2)+x(n-l))/2>  ctc'  Thc  Bickc  1 -Hodges 
estimator,  the  median  of  these  numbers,  does  not  have  the  CMVP 
(Magel  (1982)).  It  is  also  shown  there  that  the  multiply 
folded  medians  (cf.  Andrews  et  al.  (1972)  for  a  description) 
are  not  CMs. 

Magel  (1982)  has  also  shown  that  the  skipped  means  and 
the  one  step  Hubers  do  not  possess  the  CMVP.  (Andrews  et  al. 
(1972)  also  discuss  these  estimators.) 

Since  so  few  robust  estimators  of  location  have  the  CMVP, 
we  consider  algorithms  that  do  not  require  this  property. 

Leurgans  (1982)  considered  a  SLOCOM  algorithm  with  linear 
functions  of  order  statistics,  but  noted  that  the  resulting 


estimators  need  not  he  nondecreasing.  Robertson  and  Wright 
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(1974,  Example  5)  considered  the  estimation  of  ordered  modes  and 

observed  that  the  max-min  formula  applied  to  consistent 

estimators  which  are  not  CMs  may  produce  estimator  which  are 

not  consistent.  That  example  also  shows  that  the  lower 

sets  algorithm  (cf.  Barlow  et  al.  (1972))  has  this  same 

difficulty.  The  problem  stems  from  the  fact  that,  with 

estimators  that  are  not  CMs,  the  initial  estimates  may  be 

nondecreasing,  but  yet  these  algorithms  may  moditv  them 

substantially.  This  does  not  happen  with  the  I’AVA.  In 

fact,  it  is  clear  that  if  Oj  <  0,  ^  ...  <  0. ,  if  the 

estimator  T  is  consistent  and  if  min,  ,  n.  -  a.  then  0. 

l_<i<k  i  ’  i 

is  consistent  for  6^  for  i=l,2,...,k.  So  we  recommend 
the  PAY A  when  computing  order  restricted  estimates  based 
on  initial  estimators  which  are  not  CMs.  However,  there  is 
one  difficulty  with  this  approach.  The  estimates  may 
depend  on  the  order  in  which  violators  are  pooled.  In 
the  Monte  Carlo  study  that  was  conducted,  pooling  was 
always  from  left  to  right,  starting  with  the  first  violators, 
the  ith  and  i+lst  samples  are  pooled,  then 
(niXi  +  ni+l^i+1)/(ni  +  nj+i^  is  compared  with  \j  +  ,,  etc. 

3.  MONTH  CARLO  RESULTS.  Since  the  small  sample  properties 
of  order  restricted  estimators  have  proved  to  be  quite 
intractable,  even  in  the  case  of  normal  means,  a  Monte 
Carlo  study  was  conducted  to  assess  the  per forinancc  of  such 
estimators  for  various  choices  of  T.  As  was  mentioned  in 
the  last  section,  even  for  moderate  values  of  k,  several 


1  (1 

values  of  T  must  be  computed  when  employing  the  1’AVA  and  so  we 

have  not  considered  some  estimators  because  of  computational 

complexities.  IVe  have  studied  the  following  T:  the  mean; 

median;  Tukey's  trimean;  Gastwirth's  estimator  least ! ; 

symmetrically  trimmed  means,  trimming  15#  and  2 .i (  on  each 

side  (Trim  15,  Trim  25);  and  Huber's  M-estimator  with 

c  =  1.5  and  2.0  (11-1.5,  11-2.0).  IV e  have  assumed  that  t  he 

\. .  -  0.  arc  i id  with  the  distributions  studied  in  Andrews 
i  J  i 

et  al.  11972).  A  key  to  the  distributions  used  is  given  in 
Table  1.  Many  of  these  distributions  are  mixtures  and 
N ( 0 , 1 ) /U ( 0 , 1 )  represents  the  distribution  of  the  ratio  of 
independent  variables,  one  a  standard  normal  and  the  other 
a  uniform  variable  on  (0,1).  l-'or  k  =  2,5,5  and  various 
choices  of  0^  <  07  <  ...  £  0^ ;  n^  =  n,  =  ...  =  n^  =  n  =  lu 
or  20;  each  distribution  in  Table  1;  and  each  T  given  above, 

the  mean  square  error  (MSli)  associated  with  estimating  •  .  , 

__  ■> 

ie.  E(0^-9^)“,  was  estimated  based  on  5000  iterations.  The 
total  MSE,  j  .  _  ^E  ( 0  ^  -  0  )  "  ,  was  estimated  by  summing  these  K 
values  and  with  k,0^  <  0,  <  ...  0^,  and  an  error  distri¬ 

bution  fixed,  an  estimated  relative  efficiency  for  a  particular 
T  was  computed  by  taking  the  reciprocal  of  the  ratio  of  its 
total  MSE  to  the  smallest  total  MSE  for  all  T  in  the  study. 

These  estimated  relative  efficiencies  are  given  in  Tables 
II  through  V. 

In  Tables  II  and  IV,  with  k  and  n  fixed,  the  dispersions 
in  the  0  vector,  0^  -  0^,  are  varied.  However  it  is  interesting 


to  note  that  there  are  no  unusual  changes  in  the  efficiencies 
This  would  lead  us  to  believe  that  the  conclusions  that 
follow  are  valid  for  a  reasonable  range  of  dispersions  in  v . 
(Since  the  estimators,  T,  considered  are  location  statistics, 

the  0.  -  o.  have  distributions  that  are  invariant  undei 

1  l 

shifts  by  a  constant  vector.) 

First,  we  consider  the  estimators  in  groups.  I  he  median 
trimean  and  (last  are  easily  computed  h- stat  i  st  i  cs  .  (.ast 
dominates  the  median  for  all  but  the  very  heavy  tailed 
distributions  (11  and  in  some  cases  7).  However,  the  loss 
in  efficiency  is  at  most  10  percentage  points  in  all  the 
cases  considered  and  the  gain  is  as  much  as  24  percentage 
points  in  some  light  tailed  cases.  The  trimean  is  on  the 
other  extreme,  performing  even  better  for  light  tailed 
distributions  (a  gain  of  5  percentage  points  or  so  in 
efficiency),  but  considerably  worse  for  heavy  tailed  distii 
butions,  with  a  loss  of  approximately  20  percentage  points. 

In  this  group  we  recommend  Gastwirth’s  estimator  unless  the 
user  is  quite  certain  that  very  heavy  tailed  errors  are  not 
present,  then  the  trimean  could  be  used. 

Andrews  et  al.  (1972)  suggest  that  Hubers,  with  e  in  the 
range  we  have  considered,  might  be  used  in  practical 
situations.  (They  consider  the  Cauchy  distribution  to  he 
unreasonable.)  In  the  cases  considered  here,  the  efficiency 
of  H-2.0  for  a  normal  model  was  about  3®  more  than  11-1.5, 
but  that  is  reversed  in  even  the  lightest  contamination 
considered  and  there  is  as  much  as  10  percentage  points 


difference  for  moderate  tailed  errors.  So  we  recommend  ii  1  .  5 
in  this  family.  If  very  heavy  tailed  distributions  ('',11  i  can 
he  ruled  out,  an  experimenter  may  prefer  li-1.5  over  the 
other  T  in  this  study.  It  holds  it  efficiency  in  the  normal 
model,  performs  well  for  all  hut  the  very  heavy  tailed 
distributions  and  lias  the  advantages  of  a  CM. 

Among  the  trimmed  means  we  recommend  a  trimming 
proportion  around  25"  on  each  side.  The  Trim  15  has 
percentage  points  more  efficiency  for  the  normal  model, 
but  the  losses  for  moderate  and  heavy  tailed  distributions 
(4-5  and  8-10)  can  he  far  greater. 

If  one  wants  an  estimator  that  performs  well  over  tin. 
range  of  situations  considered  here,  then  the  Trim  25  and 
the  Cast  should  he  considered.  The  efficiency  of  the  (last 
for  a  normal  model  is  about  4  percentage  points  smaller, 
hut  is  up  to  7  percentage  points  larger  for  the  Cauchy 
distribution  (cf.  k  =  5) .  Castwirth's  estimator  is  easiet 
to  compute.  If  one  feels  it  is  not  necessary  to  guard 
against  very  heavy  tailed  distributions  (7,11),  then 
11-1.5  is  recommended.  It  is  about  95"  efficient  in  the 
normal  case  and  has  the  advantage  of  a  CM. 

While  it  is  not  possible  to  consider  enough  k  and  o 
to  draw  definitive  conclusions,  the  fact  tha*  the 
recommendations  given  above  arc  supported  by  each  case 
considered  strengthens  their  credibility.  Also  for  each  k 
and  0  considered,  the  Monte  Carlo  experiment  was  carried 
out  for  n  =  10  and  20, and  for  k  =  5  additional  u,  such 


as  (.0,0, 0,0,1]  and  (2,4,8,10,32),  were  considered, 
case  substantiated  the  recommendations  above. 


I.ac  h 
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TABLE  I.  KEY  TO  THE  ERROR  DISTRIBUTION  USED 


1 . 

N(0 

,D 

2  _ 

9  0% 

N  ( 0 

,D 

and 

m 

N  ( 0 , 9 1 

3 . 

Double 

Exponent ial 

4. 

50°o 

N  CO 

,D 

and 

50  o 

N(0,9j 

5. 

2  5  o 

N  C  0 

,D 

and 

7  5° 

N  (.  0 , 9  ) 

6. 

90  % 

N  (0 

,D 

and 

10% 

N(0, 100) 

7. 

7  51 

N  ( 0 

,D 

and 

25% 

N  (0,100) 

8. 

901 

N  (0 

,D 

and 

10% 

N(0,1)/U((I, 

1  j 

9. 

7  5  % 

N  (0 

,D 

and 

251 

N(0,1)/U(0, 

1) 

10. 

90% 

N  (0 

,D 

and 

10% 

X(0,1)/U(0, 

1 1 .  Cauchy 
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TABLE  II:  RELATIVE  EFFICIENCIES:  k=2,  n=20 

0  =  10, Oj 


IH  s  t  r i hut  ion 

Mean 

Med ian 

Trimean 

Cast 

11-1.5 

11-2.0 

T  r  i  m  1  5 

1  i  i  in 

1 

1 . 000 

.675 

.  86  5 

.819 

.  959 

.  987 

.916 

.  84  5 

y 

.  7  30 

.  747 

.  944 

.899 

1  .  000 

.  96  9 

.  98  9 

.  9  3  1 

3 

.  503 

.  639 

.  954 

.  883 

.  953 

.  825 

1  .  (1(10 

.925 

4 

.  642 

.  914 

.  928 

.  986 

.  856 

.773 

.  901 

1  .  Olid 

5 

.  871 

.  937 

.  955 

.  987 

.  929 

.  904 

.  94  8 

1  .  000 

6 

.  136 

.  827 

.  996 

.97  9 

.  984 

.  882 

.  984 

1  .  000 

7 

.100 

.  968 

.  724 

1 . 000 

.  763 

.  584 

.  566 

.  98  4 

8 

<.  1 

.770 

.950 

.  906 

1 . 000 

.  973 

.  98  8 

.  9  3  8 

9 

<.  1 

.828 

.987 

.965 

.  982 

.  900 

1  .  0  00 

.  9  9  1 

10 

<•1 

.  808 

.  987 

.955 

.991 

.  904 

1 . 0  00 

.  984 

1  l 

<  .  1 

1 . 000 

.  696 

.904 

.64  7 

.52  5 

.59  1 

.892 

e  = 

(  -  .  3  ,  . 

3) 

1 

1.000 

.7  03 

.882 

.  838 

.  960 

.  98  8 

.  9  2  6 

.  8  59 

2 

.751 

.769 

.953 

.  914 

1 . 000 

.971 

.  991 

.  9  3  9 

3 

.  528 

.683 

.964 

.927 

.  967 

.  860 

1.000 

.  9  5  3 

4 

.671 

.  931 

.936 

.987 

.  871 

.  795 

.  909 

1  .  00  0 

5 

.  884 

.  94  8 

.  963 

.989 

.  936 

.  914 

.  954 

1  .  000 

6 

.158 

.845 

.  988 

.982 

.977 

.882 

.  977 

1  .  OOd 

7 

.116 

.  982 

.744 

1.000 

.  785 

.614 

.  589 

.  98  2 

8 

<.l 

.  782 

.951 

.914 

1.000 

.  974 

.  988 

.93  9 

9 

<.  1 

.845 

.991 

.970 

.  987 

.  911 

1 . 000 

.  990 

10 

<.  1 

.  839 

.  993 

.  971 

.  994 

.  91  8 

1 . 000 

.  9  9  1 

11 

<.  1 

1 . 000 

.715 

.  901 

.6  73 

.5  56 

.  6  1  5 

.  8  8!) 

r- 

i  (, 


tabu: 

11:  Rlil.ATIVl:  Id 

MClIiNC 

ICS:  k= 

3,  e=(o 

,‘M  ) 

n 

=  10 

H  i  s  t  r  i  Ini  t 

ion  Mean 

Med  i  an 

T r  i  mean 

Cast 

11  1  .  5 

11-2.0 

Trim  15 

1  r  i  in  2 

1 

1  .  000 

.7  29 

.  892 

.  833 

.  9  54 

.  981 

.  9  3  9 

.  87  0 

i 

.  762 

.810 

.  955 

.  921 

1  .  000 

.  978 

.  98  5 

.  8  9  9 

3 

.  539 

.  800 

.  974 

.  96  2 

1  .  000 

.890 

.941 

.  9  7  0 

4 

.  "23 

.954 

.  94  5 

1 . 000 

.917 

.851 

.  92  0 

.  9  9  2 

3 

.8  93 

.  94  2 

.  958 

.985 

.  94  9 

.928 

.  9  7  5 

1  .  000 

6 

.173 

.  893 

.931 

.  984 

1  .  000 

.978 

.  8  0  9 

.981 

“T 

.128 

1  .  000 

.  576 

.915 

.  767 

.610 

.  4  36 

.  "64 

8 

<  .  1 

.811 

.  956 

.912 

1 . 000 

.  986 

.95“ 

.  9  3  S 

9 

•  .  1 

.  88  3 

.  96  8 

.976 

1  .  000 

.  935 

.  84  2 

.  9  9  “ 

1U 

■  .  1 

.871 

.954 

.  96  8 

1  .  000 

.  936 

.  84  4 

.97  1 

1  1 

-  .  I 

1  .  000 

.  6  3  8 

.  904 

.713 

.  5  9  7 

.  4  2  5 

.825 

n 

=  20 

1 

1 . 000 

.682 

.87  3 

.  820 

.  96  1 

.991 

.  9  3  2 

.  860 

■> 

.  7  36 

.  761 

.  949 

.  9  1  0 

1  .  000 

.  9  ~  1 

.  9  9  1 

.  9  3  8 

3 

.  511 

.  66  7 

.  978 

.914 

.997 

.8  58 

I  .  000 

.  93  4 

4 

.667 

.914 

.  923 

.  98  2 

.  863 

.“87 

.  904 

1  .  000 

5 

.837 

.919 

.92  3 

.  94  9 

.  894 

.  86  8 

.  950 

1  .  000 

0 

.  163 

.  841 

.99  2 

.  97  8 

1  .  000 

.  904 

.  984 

.  9  9  “ 

7 

.  121 

.  961 

.772 

1 . 00  0 

.  8  06 

.  6  3  3 

.585 

.  96  3 

8 

<.  1 

.  760 

.  944 

.  895 

1 . 000 

.  98  1 

.  962 

.  9  1  3 

9 

<  .  1 

.783 

.  929 

.  906 

.  9  3  9 

.  8  “  0 

1  .  000 

.  98  8 

1U 

<  .  1 

.  808 

.974 

.  96  0 

1 . 000 

.927 

.  996 

.  9  1  5 

11 

.  1 

1  .  000 

.7  24 

.  9  1  9 

.687 

.  5  "  2 

.  6  1  7 

.  8  86 

/ 
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TABU!  IV:  Rlil.AT  I VI:  FFMCI  li.NIC  1 1!S  ,  k=3,  n=20 

0  =  (-2,0,2) 


Pis  trihut  ion 

Mean 

Median 

Tr imcan 

Cast 

li-  1  .  5 

II 

-2.0 

Trim  15 

1  r  i  m  2 

l 

1 . 000 

.  681 

.87  3 

.  824 

.  96  2 

.  990 

.928 

.8  58 

2 

.734 

.  762 

.951 

.913 

1 . 000 

.  97  0 

.  9  9  1 

.  9  3  9 

3 

.  513 

.681 

.  979 

.932 

1 . 000 

.865 

.  987 

.  9  4  0 

4 

.  64  7 

.  9  0  5 

.  909 

.  96  5 

.  8  5  0 

.772 

.  8  9  9 

1  .  000 

5 

.855 

.  934 

.  950 

.971 

.915 

.  888 

.  9  5  2 

1  .  000 

(i 

.  147 

.  84  5 

.  984 

.  981 

.  99  3 

.  8  9  5 

.  96  3 

1  .  000 

- 

.116 

.  989 

.  ^06 

1.000 

.  807 

.62  2 

.529 

.  9  5  ~ 

8 

<.  1 

.7  59 

.  94  3 

.  899 

1  .  000 

.  982 

.  972 

.022 

9 

<.  1 

.  804 

.954 

.929 

.  963 

.  8  9  1 

1  .  000 

.  9  9  1 

10 

<.  1 

.818 

.  98  0 

.  959 

1  .  000 

.9  24 

.  986 

.  9  ~i  8 

LI 

-• .  1 

L .  000 

.  694 

.  903 

.680 

.  56  0 

.  590 

.888 

0  = 

( - . 33  , 

-  .  09  ,  .  3 

3) 

1 

1 . 000 

.714 

.885 

.  84  0 

.  963 

.988 

.  93  7 

.  8"2 

2 

.  762 

.  785 

.95  5 

.919 

1  .  000 

.  97  4 

.  9  9  1 

.  94  4 

3 

.555 

.699 

.  980 

.932 

1 . 000 

.  875 

1 . 000 

.  94  4 

4 

.67  7 

.926 

.932 

.  989 

.  871 

.  7  96 

.  908 

1  .  000 

5 

.825 

.911 

.913 

.  94  3 

.883 

.  856 

.  94  8 

1  .  000 

6 

.  1"1 

.  858 

.  995 

.984 

1  .000 

.  913 

.  981 

.  98  7 

7 

.118 

.  963 

.  795 

1.000 

.815 

.  646 

.602 

.  9  7  0 

8 

<•  1 

.  785 

.  950 

.908 

1.000 

.  927 

.  955 

.912 

9 

<.  1 

.  808 

.  944 

.92  2 

.  955 

.  889 

1  .  0  0  0 

.98  9 

10 

<.  1 

.  834 

.983 

.  967 

1  .  000 

.9  33 

1  .  000 

.  9  90 

11 

<.  t 

1  .  000 

.  728 

.923 

.688 

.  5  74 

.  (i  1  5 

.  8  "6 
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TABUi  V: 

RliLAT  1  VI: 

LPFICIENCIHS, 

k  =  5,  n= 

10 

0  = 

(-.8,-. 

3 , 0  ,  .  9  , 

1.5) 

t  ribut ion 

Mean 

Median 

Trimean 

(.as  t. 

11-  1  .  5 

H-2.0 

T  rim  1 

1 

1 .  000 

.757 

.  903 

.  849 

.  96  2 

.  985 

.  9  4  9 

-> 

.782 

.828 

.  958 

.  92  4 

1  .  000 

.  982 

.  987 

3 

.  512 

.  668 

.979 

.916 

1  .  000 

.  86  1 

.  95  4 

4 

.725 

.  96  2 

.  9.36 

.991 

.911 

.  8  4  (> 

.  9  2  5 

5 

.  884 

.  94  5 

.951 

.  97  3 

.938 

.  9  1  () 

.  97  1 

0 

.  197 

.  901 

.  95  1 

.  978 

1 . 000 

.  926 

.  8  3  1 

7 

.  145 

1 . 000 

.645 

.  94  5 

.823 

.67  2 

.  509 

8 

<.  1 

.826 

.957 

.921 

1.000 

.988 

.  961 

9 

<•  1 

.888 

.  967 

.  971 

.  987 

.929 

.  92  0 

10 

<.  1 

.878 

.  948 

.  961 

1 . 000 

.  945 

.  886 

11 

<.  1 

1.000 

.724 

.919 

.687 

.  572 

.430 

0  = 

(-2.2,- 

.7,0,0, 

•  7) 

1 

1.000 

.  742 

.8  96 

.  839 

.  958 

.  985 

.  938 

2 

.  776 

.826 

.  955 

.  925 

1 . 000 

.  98  3 

.  98  2 

3 

.  54  3 

.796 

.954 

.938 

1  .  000 

.891 

.  9  1  0 

4 

.723 

.  95  9 

.935 

.  990 

.  912 

.  846 

.  92  1 

5 

.889 

.  947 

.956 

.980 

.94  2 

.921 

.  9  7  5 

6 

.192 

.  888 

.942 

.  973 

1  .  000 

.  924 

.801 

7 

.  147 

1.000 

.6  30 

.936 

.  823 

.670 

.  500 

8 

<.  1 

.820 

.955 

.917 

1 . 000 

.  989 

.  95  7 

9 

<.  1 

.  878 

.957 

.  960 

.  978 

.920 

.912 

10 

<.  1 

.871 

.943 

.958 

1 . 000 

.  941 

.857 

11 

<.  1 

1.000 

.642 

.  912 

.743 

.632 

.431 

T rim  2  5 
.889 
.  9  5  4 
.  990 
1  .000 
1  .  000 
.  975 
.8  54 
.  94  0 
1  .  000 
.  <48  7 
.  8  50 

.  875 
.  953 
.942 
1  .  000 
1  .  000 
.  9  5  8 
.840 
.  94  3 
1  .  000 
.  90  9 
.  834 


